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Abstract. In view of the recent results from the MiniBooNE experiment we revisit the global 
neutrino oscillation fit to short-baseline neutrino data by adding one, two or three sterile 
neutrinos with eV-scale masses to the three Standard Model neutrinos [1]. We find that four- 
neutrino oscillations of the (3+1) type, which have been only marginally allowed before the 
recent MiniBooNE results, become even more disfavored with the new data. In the framework 
of so-called (3+2) five-neutrino mass schemes the MiniBooNE results can be nicely reconciled 
with the LSND appearance evidence thanks to the possibility of CP violation available in such 
oscillation schemes; however, the tension between appearance and disappearance experiments 
represents a serious problem in (3+2) schemes, so that these models are ultimately not viable. 
This tension remains also when a third sterile neutrino is added, and we do not find a significant 
improvement of the global fit in a (3+3) scheme. 



1. Introduction 

Recently the first results from the MiniBooNE (MB) experiment [2, 3] at Fermilab have been 
released on a search for — > v e appearance with a baseline of 540 m and a mean neutrino energy 
of about 700 MeV. The primary purpose of this experiment is to test the evidence of Un — * v e 
transitions reported by the LSND experiment at Los Alamos [4] with a very similar L/E range. 
Reconciling the LSND signal with the other evidence for neutrino oscillations is a long-standing 
challenge for neutrino phenomenology, since the mass-squared differences required to explain the 
solar, atmospheric and LSND experimental results in terms of neutrino oscillations differ from 
one another by various orders of magnitude. Consequently, there is no consistent way to explain 
all these three signals invoking only oscillations among the three known neutrinos. Therefore, in 
order to explain the LSND anomaly one had to invoke an extension of the three-neutrino mixing 
scenario, introducing either a mechanism to generate at least a third mass-square difference, 
or a new form of flavor transition beyond oscillations. Following Ref. [1], in this talk we will 
concentrate on the first possibility, starting from models with one extra sterile neutrino (Sec. E]) 
and then considering models with two and three sterile neutrino states (Sec. [3]). 

2. Four-neutrino mixing 

In four-neutrino models, one extra sterile state is added to the three weakly interacting ones. 
The relation between the flavor and the mass eigenstates can be described in terms of a 4 x 4 
unitary matrix U, which generalizes the usual 3x3 leptonic matrix of the Standard Model. 
There are six possible four-neutrino schemes that can accommodate the results from solar and 
atmospheric neutrino experiments and contain a third much larger Am 2 . They can be divided 




into two classes: (3+1) and (2+2). In the (3+1) schemes, there is a group of three close-by 
neutrino masses that is separated from the fourth one by the larger gap. In (2+2) schemes, 
there are two pairs of close masses separated by the large gap. While different schemes within 
the same class are presently indistinguishable, schemes belonging to different classes lead to very 
different phenomenological scenarios. 

A characteristic feature of (2+2) schemes is that the extra sterile state cannot be 
simultaneously decoupled from both solar and atmospheric oscillations. To understand why, 
let us define rj s = J2iesoL \ U s i\ 2 and c s = X^'gatm l^sjl 2 ' where the sums in i and j run over 
mass eigenstates involved in solar and atmospheric neutrino oscillations, respectively. Clearly, 
the quantities rj s and c s describe the fraction of sterile neutrino relevant for each class of 
experiment. Results from atmospheric and solar neutrino data imply that in both kind of 
experiments oscillation takes place mainly between active neutrinos. Specifically, from Fig. 46 
of Ref. [5] we get r/ s < 0.31 and c s < 0.36 at the 3<r level. However, in (2+2) schemes unitarity 
implies rj s + c s = 1. A statistical analysis using the parameter goodness of fit (PG) proposed 
in [6] gives Xpg = 30.7 for 1 d.o.f., corresponding to a 5.5a rejection (PG = 3 x 10~ 8 ) of the 
(2+2) hypothesis. These models are therefore ruled out at a very high confidence level, and in 
the rest of this talk we will not consider them anymore. 

On the other hand, (3+1) schemes are not affected by this problem. Although the 
experimental bounds on rj s and c s quoted above still hold, the condition rj s + c s = 1 no longer 
applies. For what concerns neutrino oscillations, in (3+1) models the mixing between the sterile 
neutrino and the three active ones can be reduced at will, and in particular it is possible to 
recover the usual three-neutrino scenario as a limiting case. However, as widely discussed in 
the literature (see, e.g., Ref. [7] and references therein) these models are strongly disfavored as 
an explanation of LSND by the data from other short-baseline (SBL) laboratory experiments. 
In the limit A?n,L SND ^> Am| TM 3> Amg 0L the probability P v -> Ve which is relevant for LSND 
as well as for KARMEN [8], NOMAD [9] and MiniBooNE is driven by the large Am^, and is 
given by 

where L is the distance between source and detector. The LSND, KARMEN, NOMAD 
and MiniBooNE experiments give allowed regions in the (Am| 1; {Ue^U^l 2 ) plane which can 



be directly obtained from the corresponding two-neutrino exclusion plots [3,4,8,9]. At the 
light of the recent MiniBooNE result which is consistent with no oscillations above 475 MeV, 
practically all the LSND region is now excluded. In addition, further constraints on It/^t/^l 2 
can be obtained by combining together the bounds on |i/ e 4| and |c7^ 4 | derived from reactor and 
accelerator experiments (mainly Bugey [10] and CDHS [11]) as well as solar and atmospheric 
data. Present bounds are plotted in the two leftmost panels of Fig. [T] as a function of Am^; 
for a discussion on the sensitivity of future experiments, see Ref. [12]. The results of the global 
analysis presented in Ref. [1], which includes atmospheric and long-baseline data together with 
all the short-baseline experiments observing no evidence (NEV), are summarized in the right 
panel of Fig. [TJ Using the PG test discussed above we find Xpg = 24.7 for 2 d.o.f., corresponding 
to a 4.6o" rejection (PG = 4 x 10~ 6 ) of the (3+1) hypothesis. These results show that (3+1) 
schemes are now ruled out as a possible explanation of LSND [1]. In addition, it should be noted 
that the low-energy excess observed by MiniBooNE at E u < 475 MeV cannot be explained in 
terms of oscillations with only one large mass-squared difference, thus adding another problem 
to these models in case this excess is confirmed to be a real signal. 

3. Five-neutrino and six-neutrino mixing 

Five-neutrino schemes of the (3+2) type are a straight-forward extension of (3+1) schemes. In 
addition to the cluster of the three neutrino mass states accounting for "solar" and "atmospheric" 
mass splittings now two states at the eV scale are added, with a small admixture of v e and 
to account for the LSND signal. In the Appendix of Ref. [13] it was suggested that such models 
could somewhat relax the tension existing between short-baseline experiments and the LSND 
data. In Ref. [14] a complete analysis was performed, finding that indeed the disagreement 
between LSND and null-result experiments is reduced. Here we will reconsider this possibility 
at the light of the new MiniBooNE data. As explained in Ref. [3] , MiniBooNE found no evidence 
of oscillations above 475 MeV, whereas below this energy a 3.6<r excess of 96 ± 17 ± 20 events is 
observed. Whether this excess comes indeed from i/„ — > v e transitions or has some other origin 
is under investigation [3]. Lacking any explanation in terms of backgrounds or systematical 
uncertainties, we will present the results obtained using both the full energy range from 300 MeV 
to 3 GeV ("MB300") and for the restricted range from 475 MeV to 3 GeV ("MB475"). 

As for (3+1) models, in (3+2) schemes the appearance data (LSND, KARMEN, NOMAD, 
and MiniBooNE) can be described using the SBL approximation Amg 0L pa and Am 2 TM pa 0, 
in which case the relevant transition probability is given by 

Pvn-tve = 4 \Ue4U^\ 2 SIV? 41 + 4 \U e5 U^ 5 \ 2 sin 2 51 

+ 8|£/ e4 t/ At 4| {UesU^ I sin 04i sin 51 cos (054 - 5) , (2) 

with the definitions (pij = /^m 2 -L/AE and 5 = arg([/* 4 L^ 4 [/e5[/* 5 ). Eq. ([2]) holds for neutrinos 
(NOMAD and MB); for anti-neutrinos (LSND and KARMEN) one has to replace 5 -5. 
Note that Eq. ([2]) is invariant under the transformation 4^5 and 5 <-> —5, and depends 
only on the combinations {Ue^U^l and \U e ^U^\. An important observation is that non-trivial 
values of the complex phase 5 lead to CP violation, and hence in (3+2) schemes much more 
flexibility is available to accommodate the results of LSND (anti- neutrinos) and MB (neutrinos). 
In Fig. [2] we show the prediction for MB and LSND at the best fit points of the combined MB, 
LSND, KARMEN, NOMAD analysis. As can be seen from this figure, MB data can be fitted 
very well while simultaneously explaining the LSND evidence. Furthermore, in this case also 
the low energy MB data can be explained, and therefore, in contrast to (3+1) schemes, (3+2) 
oscillations offer an appealing possibility to account for this excess. For both MB475 and MB300 
a goodness-of-fit of 85% is obtained, showing that MB is in very good agreement with global 
SBL appearance data including LSND. 
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Figure 2. Spectral data for the MiniBooNE (left) and LSND (right) experiments, calculated at 
the best fit point of a combined fit of LSND, KARMEN, NOMAD and MB in (3+2) schemes. 



On the other hand, once disappearance data are included in the analysis, the quality of the fit 
decreases considerably. Indeed, even in (3+2) schemes short-baseline experiments pose stringent 
bounds on the mixing angles \U e i\ and \Uai\, in close analogy with (3+1) models described in 
Sec. Eland shown in the two leftmost panels of Fig. [TJ Since rather large values of lUe^U^l and 
{UesUusl are needed to account for the negative result of MiniBooNE as well as the positive signal 
of LSND, one expects that reconciling appearance and disappearance data will be a problem also 
within (3+2) models. This tension is illustrated in Fig. where the projections of the allowed 
regions in the plane of the appearance amplitudes \U e iU^\ and \U e ^U^\ are shown. Indeed the 
opposite trend of the two data sets is clearly visible, especially when the low energy excess in 
MB is included (right panel). In order to quantify this disagreement one can apply the PG test 
to appearance versus disappearance data without MB, with MB475, and with MB300: 



APP vs DIS: ^ 



Xp G = 17.5 , PG = 1.5 x l(T d (no MB), 

Xp G = 17.2 , PG = 1.8 x 1(T 3 (MB475), (3) 

Xp G = 25.1, PG = 4.8xlO" 5 (MB300). 



From these numbers we conclude that also in (3+2) schemes the tension between appearance 
and disappearance experiments is quite severe. If MB475 is used the result is very similar to 
the pre-MiniBooNE situation implying inconsistency at about 3. la, whereas in case of the full 
MB300 data the tension becomes significantly worse (about 4a), since appearance data are 
more constraining because of the need to accommodate LSND as well as the MB excess at low 
energies. 

Finally, since there are three active neutrinos it seems natural to consider also the case of 
three sterile neutrinos. If all three additional neutrino states have masses in the eV range and 
mixings as relevant for the SBL experiments under consideration, such a model will certainly 
have severe difficulties to accommodate standard cosmology [15]. Besides this fact, the results 
of the search performed in Ref. [1] show that there is only a marginal improvement of the fit by 
1.7 units in \ 2 f° r MB475 (3.5 units for MB300) with respect to (3+2), to be compared with 
four additional parameters in the model. Hence, the conclusion is that here are no qualitatively 
new effects in the (3+3) scheme. The conflict between appearance and disappearance data 
remains a problem, and the additional freedom introduced by the new parameters does not 
relax significantly this tension. 
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Figure 3. Allowed regions at 90% and 99% CL in (3+2) schemes from the analysis of appearance 
and disappearance data, projected onto the plane of [U^U^l and \U e ^,U^\. 



4. Conclusions 

We have considered the global fit to short-baseline neutrino oscillation data including the recent 
data from MiniBooNE, in the framework of (3+1), (3+2) and (3+3) oscillation models. Four- 
neutrino models are ruled out since (a) the don't allow to account for the low energy event excess 
in MB, (b) MiniBooNE result cannot be reconciled with LSND, and (c) there is severe tension 
between appearance and disappearance experiments. Five-neutrino models provide a nice way 
out for problems (a) and (b), but fail to resolve (c). Similarly, six-neutrino models do not offer 
qualitatively new effects with respect to (3+2). In all cases we find severe tension between 
different sub-samples of the data, hence we conclude that at the light of present experimental 
results it is not possible to explain the LSND evidence in terms of sterile neutrinos. 
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